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Abstract 

The statistical properties of trajectories of eigenvalues of Gaussian 
complex matrices whose Hermitean condition is progressively broken 
are investigated. It is shown how the ordering on the real axis of the 
real eigenvalues is reflected in the structure of the trajectories and also 
in the final distribution of the eigenvalues in the complex plane. 

1 Introduction 

A few years after the introduction by E. Wigner[l] of the Gaussian ensemble 
of Hermitean random matrices as a basis of a statistical theory of spectra of 
complex many-body systems, Ginibre investigated properties of the Gaussian 
matrices when the Hermitian condition is removed[2j. The Hermitian ones 
have found applications [3], especially since the link with the characterization 
of the manifestation of chaos in quantum mechanics has been established |4]. 
More recently, the nonHermitian class has attracted great interest [5]. 
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The interpolation between universality classes of matrices is a standard 
subject of investigation in random matrix tlieories|S]. The Wigner-Ginibre 
transition, which is our present object of study, has already been matter 
of investigation. It has been established that asymptotically, its density of 
eigenvalues, starting from Wigner's semi-circle law distribution on the real 
axis, evolves as the transition proceeds through an elliptic shape, ending up 
in an uniform circular distribution at the Ginibre situation [7]. In addition, 
an asymptotic weak non-Hermiticity regime [8] with important applications 
to quantum open systems [9] has been found in the intermediate regime of 
transition. 

From the physical point of view, this transition can be seen as the passage 
of a gas from a ID configuration to an isotropic 2D configuration. In fact, it 
is well known that the positions of eigenvalues of Gaussian ensembles can be 
considered as those of charges of a Coulomb gas. It also has been recognized 
that the real eigenvalues of the Hermitian matrices have properties of a ID 
gas of bosons, the so-called Girardeau gas[10]. On the other hand, the joint 
distribution of the complex eigenvalues of the complete non-Hermitian ma- 
trices, that is, the Ginibre matrices have the same structure of the Laughlin 
wave function of a 2D gas|llj. 

The behavior of individual eigenvalues of Hermitian random matrices has 
been a matter of investigation since the pioneering works which described 
the distribution of largest eigenvalues, an achievement followed by several 
applications [12]. Recently, the order statistics problem of obtaining the dis- 
tribution of all eigenvalues considered as an ordered sequence of random 
variables has been addressed [13]. Here we are interested in studying how 
that ordering on the real axis reflects on the eigenvalue trajectories in the 
complex plane as the Hermitian condition is progressively removed and also 
in the final eigenvalue distribution. In other words, do the complex eigen- 
values retain a memory of their initial positions on the real axis? To do this 
investigation, we resort to a system of differential equations which describes 
the motion of the eigenvalues as a function of the parameter breaking the 
Hermitian condition. 

2 Transition equations and results 

The Ginibre ensemble consists of complex random matrices S whose joint 
distribution of elements is given by 
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P{S) = exp [-tr(>S^5)] (1) 

in which no Hermitian condition is imposed. Speciahzing to the case we are 
interested, namely complex matrices of size A^, their eigenvalues, for large 
values of A^, are uniformly distributed in a disk of radius Beyond this 
circle, at a radius r, the density decays as[ 



2u 

where u = r — \/iV, such that eigenvalues can be found up to ~ \/N + 2. 

Taking a matrix out of this ensemble, we define a new matrix, H(t), by 
the relation 

where the parameter t varies from zero to one. It is relevant to mention that 
others' parametrizations have been used to study this transition [9J; however, 
for our purpose of investigating trajectories, the parameter t defined by Eq. 
([3]) is more convenient. Of course, our t is related to those others parameters 
by a simple transformation. With t > 0, Eq. ([3]) together with its adjoint, 
can be inverted to express S in terms of H and H'' as 

Substituting ([1]) in ([T]), we obtain the density distribution of the t— dependent 
matrix elements of H 



P{H) = KN{t)exp -tr 



(5) 



For if: = the matrices H{0) are Hermitian such that, in the limit t — ?■ 0, (jS]) 
becomes 

P[H{0)] = exp (-tr [H'^{0)] ) . (6) 

Therefore H{0) belongs to the Gaussian Unitary Ensemble whose eigenvalues 
are distributed on the real axis according to the Wigner semi-circle law[T3] 



3 



p{x) = ^V2N-x^. (7) 
The joint distribution of the eigenvalues apart from a constant is given by 



P(Zi, Z2, Zn) 



exp 



^ / 2\ 



n 

j>t 



Zj Zi 



P= exp [-W{zi,...,zn)] 
(8) 



where 



W{zi,...,zn) 



k=l ^ ' j>i 



(9) 



Therefore, tlie eigenvalues can be considered as positions of point charges 
under the action of a confining potential harmonic in the two Cartesian axis 
and a repulsive ID Coulomb force. 

For intermediate values of t, that is < t < 1, asymptotically, that is 
when becomes large, the eigenvalues fill an ellipse [7] whose axes for our 
parameterization are 



and 



2N 

1 + ^2 



2N 



(10) 



(11) 



1 + ^2 

As 5* = -f^(l) belongs to the Ginibre ensemble of non-Hermitian matrices, 
H{t) undergoes a transition from the Wigner to the Ginibre ensemble with 
eigenvalues moving along trajectories in the complex plane. Their motion 
is governed by a system of differential equations which can be deduced in 
the following way [15]. The matrices H{t) are diagonalized by the similarity 
transformation 



D = Q-'HQ, (12) 

where D is a diagonal matrix whose diagonal contains the complex eigenval- 
ues while Q is a matrix of size A^ which contains the eigenvectors. On the 
hand, the adjoint of Eq. (HM 
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Dt = Qt^t(g-i)t (^3) 

shows that the inverse of the adjoint of Q diagonahzes the Hermitian of H. 
All these matrices are functions of the parameter t. Taking the derivative 
(denoted by a dot) with respect to t of ( IT2!) we obtain 

b=[D,U] + P, (14) 

where 

P = Q-'HQ (15) 

and 

U = Q-'Q = -Q-'Q (16) 
while for the derivative of P we have 



P = [P, U] + Q-^HQ. (17) 
The diagonal part of (1141) gives 

Dkk = Zk = Pkk, (18) 

while the off-diagonal part 

Dki = = DkkUki - UmDu + Pki (19) 

yields 

Uki = (20) 

Zk — Zi 

Using the arbitrariness of the matrix U we impose the necessary condition 
that its diagonal elements vanish, that is Ukk = 0. Finally, from ([T6|) we 
derive the equations 

o« = Er^ (21) 

and 
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for the evolution of the eigenvectors. The above equations together with the 
initial conditions provided by the Hermitian matrix i^(0) form a complete 
system of first order differential equations which numerically can be solved 
to obtain eigenvalues and eigenvectors along the transition. 

Insofar as the dependence of H with the parameter t has not been spec- 
ified, this set of differential equations is general. Assuming this dependence 
to be given by (E]), the second derivative H vanishes due to the linear depen- 
dence with t and, explicitly, Eq. f[T7|) gives 

jS \ ^ '^PkmPmk /'OQ^ 



for the diagonal elements and 



Zk 



kl — 

m^k,mj^l 



J2 PkmPn^ki^ + ^ (24) 



for the others elements. These two above equations together with Eq. f|T9l) 
form a set of 2N[N + 1) coupled differential equations to be integrated. 
Regarding the initial conditions, since i^(0) is Hermitian, its eigenvalues are 
real and its eigenvector matrix is unitary, that is = QK Taking this into 
account, we find that initially the matrix P is given by 



such that 



whose diagonal part 



P{0) = QHO)^^Q{0), (25) 



P(0) = -pt(o), (26) 



PkkiO) = -P:M (27) 



shows that the diagonal elements of P{0) are pure imaginary. Substituting 
these initial 'velocities' in f|T8|) . we find that the eigenvalues leave the real 
axis at t = perpendicularly. 
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We observe that integrating these equations between the initial value 
t = and some final value t = tj is equivalent to diagonalize H{t) at each 
intermediate value of t. This is illustrated in Fig. 1, in which = 20 
eigenvalue trajectories obtained performing the integration from t = to t = 
1 the equations of motion are compared with the result of diagonalizing the 
matrices at several intermediate values of t. The agreement is perfect showing 
that the above system is reliable and amenable to numerical integration. 

Fig. 1 also shows that eigenvalues may have complicated trajectories in 
the complex plane. Therefore, in order to reveal structures the trajectories 
may have, it is necessary to accumulate results of many simulations, that is 
to construct an ensemble. This is done in Fig. 2 in which trajectories were 
obtained by evolving eigenvalues of 60 different initial Hermitian matrices, 
clearly exhibiting a structure. Indeed, trajectories starting on the real axis 
at the edge, cover regions of a meniscus shape whose curvature decreases as 
eigenvalues more to the center are considered. Trajectories starting at the 
central region of the initial spectrum move inside strips. The presence of 
these strips is better seen in Fig. 3 where points correspond to eigenvalues at 
the end of the transition, that is for t = 1. The structure of the distribution 
of points show that, statistically, eigenvalues preserve to some extent the 
relative positions they had on the real axis. 

3 Concluding remarks 

In conclusion, evolving under the action of the external confining harmonic 
potential and the repulsion force among them, eigenvalues present structures 
in their trajectories and in their final distribution in the complex plane which 
refiect the ordering they have on the real axis. We remark that the present 
study can be considered as an instance of the so-called parametric statistics 
used to characterize the evolution of individual eigenvalues as a function of 
an external parameter |16j. In the case of complex eigenvalues, parametric 
evolution has been experimentally investigated by considering resonances 
trapped in an open microwave cavity in which the slit width can be varied [T7]. 

This work is supported by a CAPES/COFECUB project and MPP is 
supported by the Brazilian agencies CNPq and FAPESP. 
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Figure Captions 

Fig. 1 The hnes are N = 20 eigenvalue trajectories generated integrating 
the equations of motion, and the dots over the hnes are eigenvalues obtained 
by diagonalizing the matrix. 

Fig. 2 Ensemble of trajectories of 60 matrices of size = 20 generated 
evolving eigenvalues initially at the edges, the middle parts, and the center 
of the spectra, together with the density circle of radius \/N. 

Fig. 3 Ensemble of eigenvalues at the Ginibre regime evolved from eigen- 
values at the edges, the middle parts, and the center of the spectra, together 
with the density circle of radius \/N. 
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